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ABSTRACT 

In this article, we present the results of a series of twelve 3.6-cm radio con- 
tinuum observations of T Tau Sb, one of the companions of the famous young 
stellar object T Tauri. The data were collected roughly every two months be- 
tween September 2003 and July 2005 with the Very Long Baseline Array (VLBA). 
Thanks to the remarkably accurate astrometry delivered by the VLBA, the ab- 
solute position of T Tau Sb could be measured with a precision typically better 
than about 100 micro-arcseconds at each of the twelve observed epochs. The 
trajectory of T Tau Sb on the plane of the sky could, therefore, be traced very 
precisely, and modeled as the superposition of the trigonometric parallax of the 
source and an accelerated proper motion. The best fit yields a distance to T Tau 
Sb of 147.6 ± 0.6 pc. The observed positions of T Tau Sb are in good agreement 
with recent infrared measurements, but seem to favor a somewhat longer orbital 
period than that recently reported by Duchene et al. (2006) for the T Tau Sa/T 
Tau Sb system. 
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1. Introduction 

To provide accurate observational constraints for pre-main sequence evolutionary mod- 
els, and thereby improve our understanding of star-formation, it is crucial to measure as 
accurately as possible the properties (age, mass, luminosity, etc.) of individual young stars. 
The determination of most of these parameters, however, depends critically on the often 
poorly known distance to the object under consideration. While the average distance to 
nearby low-mass star-forming regions (e.g. Taurus or p— Ophiuchus) has been estimated to 
about 20% precision using indirect methods (Elias 1978a,b, Kenyon et al. 1994, Knude & 
Hog 1998), the line-of-sight depth of these regions is largely unknown, and accurate distances 
to individual objects are still missing. Even the highly successful Hipparchos mission (Perry- 
man et al. 1997) did little to improve the situation (Bertout et al. 1999) because young stars 
are still heavily embedded in their parental clouds and are, therefore, faint in the optical 
bands observed by Hipparchos. Future space missions such as GAIA will undoubtedly have 
the capacity to accurately measure the trigonometric parallax of optically fainter stars, but 
these missions will still be unable to access very deeply embedded sources, and will only 
start to provide results in about a decade. In the meantime, extremely high quality infrared 
and X-ray surveys of many star-forming regions are being obtained (e.g. Evans et al. 2003, 
Giidel et al. 2007), and their potential cannot be fully exploited because of the unavailability 
of good distance estimates. 

Low-mass young stars often generate non-thermal continuum emission produced by the 
interaction of free electrons with the intense magnetic fields that tend to exist near their sur- 
faces (e.g. Feigelson & Montmerle 1999). Since the magnetic field strength decreases quickly 
with the distance to the stellar surface (as r~ 3 in the magnetic dipole approximation), the 
emission is strongly concentrated to the inner few stellar radii. If the magnetic field intensity 
and the electron energy are sufficient, the resulting compact radio emission can be detected 
with Very Long Baseline Interferometers (VLBI -e.g. Andre et al. 1992). The relatively 
recent possibility of accurately calibrating the phase of VLBI observations of faint, compact 
radio sources using nearby quasars makes it possible to measure the absolute position of these 
objects (or, more precisely, the angular offset between them and the calibrating quasar) to 
better than a tenth of a milli-arcsecond (Loinard et al. 2005, see also below). This level of 
precision is sufficient to constrain the trigonometric parallax of sources within a few hundred 
parsecs of the Sun (in particular of nearby young stars) with a precision better than a few 
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percents using multi-epoch VLBI observations. 

Taking advantage of this situation, we have recently initiated a large project aimed 
at accurately measuring the trigonometric parallax of a significant sample of magnetically 
active young stars in nearby star-forming regions (Taurus, p— Ophiuchus, Perseus, Serpens, 
and Cepheus) using the 10-element Very Long Baseline Array (VLBA) of the National Radio 
Astronomy Observatory (NRAO). In the present article, we will concentrate on T Tau Sb, 
one of the members of the famous young stellar system T Tauri (see e.g. Duchene et al. 2006 
for a recent summary of the properties of that system). T Tau Sb has long been known 
to be associated with a compact non-thermal radio source (Skinner & Brown 1994; Phillips 
et al. 1993, Johnston et al. 2003) characterized by strong variability and significant circular 
polarization. An extended, thermal, radio halo studied in detail by Loinard et al. (2007) and 
probably related to stellar winds, also exist around T Tau Sb. While this extended structure 
contributes to the total radio flux as measured, for instance, with the VLA, it is effectively fil- 
tered out in Very Long Baseline Interferometry experiments. Indeed, in the intercontinental 
VLBI observations published by Smith et al. (2003), only about 40% of the simultaneously 
measured VLA flux density is retrieved. The radio source detected by Smith et al. (2003) is 
very compact (R < 15 R©), and its flux was about 3 mJy at the time of their observations. 
Its trajectory over the plane of the sky was studied by Loinard et al. (2005) using a series of 
7 VLBA observations. Unfortunately, these data were recently found to have been affected 
by a bug that caused the VLBA correlator to use predicted rather than measured Earth 
Orientation Parameters (see |http:/ /www.vlba.nrao.edu/astro/messages/eop/[ ). This prob- 




lem corrupted the visibility phases, and strongly affected the quality of the astrometry of the 
data published in Loinard et al. (2005). The post-fit rms for the data published by Loinard 
et al. (2005) was about 250 mas compared with 60-90 mas for the present data (see below). 
Here, we will re-analyze these VLBA data, and combine them with 5 newer observations to 
measure the trigonometric parallax, and study the proper motion of T Tau Sb. 



In this paper, we will make use of a series of twelve continuum 3.6 cm (8.42 GHz) ob- 
servations of T Tau Sb obtained every two months between September 2003 and July 2005 
with the VLBA (Tab. 1). Our phase center was at aj 20 oo.o — 04 h 21 m 59. s 4263, <5/2ooo.o = 
+19°32'05''730, the position of the compact source detected by Smith et al. (2003). Each ob- 
servation consisted of series of cycles with two minutes spent on source, and one minute spent 
on the main phase-referencing quasar J0428+1732, located 2?6 away. J0428+1732 is a very 
compact extragalactic source whose absolute position (aj2ooo.o = 04 h 28 m 35?633679,5j2ooo.o = 
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Observations and data calibration 
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17°32'23''58799) is known to better than 1 milli-arcsecond (a a = 0.59 mas, as = 0.89 mas; 
Beasley et al. 2002). During the first 6 observations, the secondary quasar J0431+1731 was 
also observed periodically -about every 30 minutes- to check the astrometric quality of the 
data, and to compare our results with those of Smith et al. (2003) who used J0431+1731 as 
their phase calibrator. A detailed comparison with the results of Smith et al. (2003) and with 
the numerous VLA observations available from the literature, however, will be postponed to 
a forthcoming article. 

The data were edited and calibrated using the Astronomical Image Processing System 
(AIPS -Greisen 2003). The basic data reduction followed the standard VLB A procedures 
for phase- referenced observations. First, the most accurate measured Earth Orientation 
Parameters obtained from the US Naval Observatory database were applied to the data 
in order to correct the erroneous values initially used by the VLBA correlator. Second, 
dispersive delays caused by free electrons in the Earth's atmosphere were accounted for 
using estimate of the electron content of the ionosphere derived from Global Positioning 
System (GPS) measurements. A priori amplitude calibration based on the measured system 
temperatures and standard gain curves was then applied. The fourth step was to correct the 
phases for antenna parallactic angle effects, and the fifth was to remove residual instrumental 
delays caused by the VLBA electronics. This was done by measuring the delays and phase 
residuals for each antenna and IF using the fringes obtained on a strong calibrator. The final 
step of this initial calibration was to remove global frequency- and time-dependent phase 
errors using a global fringe fitting procedure on the main phase calibrator (J0428+1732), 
which was assumed at this stage to be a point source. 

In this initial calibration, the solutions from the global fringe fit were only applied to 
the main phase calibrator itself. The corresponding calibrated visibilities were then imaged, 
and several passes of self-calibration were performed to improve the overall amplitude and 
phase calibration. In the image obtained after the self-calibration iterations, the main phase 
calibrator is found to be slightly extended. To take this into account, the final global fringe 
fitting part of the reduction was repeated using the image of the main phase calibrator as 
a model instead of assuming it to be a point source. Note that a different phase calibrator 
model was produced for each epoch to account for possible small changes in the main cali- 
brator structure from epoch to epoch. The solutions obtained after repeating this final step 
were edited for bad points and applied to the target source. Using an image model for the 
calibrator rather than assuming a point source improved the position accuracy by a few tens 
of //as. 

Because of the significant overheads that were necessary to properly calibrate the data, 
only about 3 of the 6 hours of telescope time allocated to each of our observations were 
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actually spent on source. Once calibrated, the visibilities were imaged with a pixel size of 50 
/ias after weights intermediate between natural and uniform (ROBUST = in AIPS) were 
applied. This resulted in a typical r.m.s. noise level of 70 /iJy for most observations, though 
for a few epochs with less favorable weather conditions, the noise level exceeded 100 /xJy 
(Tab. 1). T Tau Sb was detected with a signal to noise better than 10 at each epoch (Tab. 
1), and its absolute position (listed in columns 2 and 4 of Tab. 1) was determined using a 
2D Gaussian fitting procedure (task JMFIT in AIPS). This task provides an estimate of the 
position error (columns 3 and 5 of Tab. 1) based on the expected theoretical astrometric 
precision of an interferometer (Condon 1997). Systematic errors, however, usually limit the 
actual precision of VLBI astrometry to several times this theoretical value (e.g. Fomalont et 
al. 1999, Pradel et al. 2006). At the frequency of the present observations, the main sources 
of systematic errors are inaccuracies in the troposphere model used, as well as clock, antenna 
and a priori source position errors. These effects combine to produce a systematic phase dif- 
ference between the calibrator and the target that limits the precision with which the target 
position can be determined. We did not attempt to correct for these systematic effects here, 
and will, therefore, assume that the true error on each measurement is the quadratic sum 
of the random error listed in Tab. 1 and a systematic contribution. The latter is difficult to 
estimate a priori, and will be deduced from the fits to the data. 



The displacement of T Tau Sb on the celestial sphere is the combination of its trigono- 
metric parallax (n) and its proper motion. For isolated sources, it is common to consider 
linear and uniform proper motions, so the right ascension (a) and the declination (5) vary 
as a function of time t as: 



where cto and Sq are the coordinates of the source at a given reference epoch, fi a and fis 
are the components of the proper motion, and f a and f$ are the projections over a and 
5, respectively, of the parallactic ellipse. The latter functions are given by (e.g. Seidelman 
1992): 



3. Astrometry fits 



a(t) 
S(t) 



«0 + (Va COS 5)t + 7lf a (t) 

S + fj,st + nfs(t), 



(1) 
(2) 
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/«(*) 
fs(t) 



(X sin ai — Y cos a±) / (15 cos Si) 
(X cos ai sin 5i + F sin cti sin #i • 



Z cos Si), 



(3) 
(4) 



where (X,Y,Z) are the barycentric coordinates of the Earth in Astronomical Units, and where 
ai = CK — 7rf a (t) and Si = S — nfsit) are the coordinates of the barycentric place of the source 
at each epoch. Note that f a and fs depend implicitly on time (through X,Y,Z) and explicitly 
on the coordinates of the source. The latter dependence on ai and Si (which are only known 
if the trigonometric parallax is known) implies that the fitting procedure must be iterative. 
The barycentric coordinates of the Earth (as well as the Julian Date of each observation) 
were calculated using the Multi-year Interactive Computer Almanac (MICA) distributed as 
a CDROM by the US Naval Observatory. They are given explicitly in Tab. 2 for all epochs. 
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Fig. 1. — Measured positions of T Tau Sb and best fit without (a) and with (b) accelera- 
tion terms. The observed positions are shown as ellipses, the size of which represents the 
magnitude of the errors. Note the very significant improvement when acceleration terms are 
included. 



- 8- 



As mentioned earlier, T Tau Sb is a member of a multiple system (e.g. Loinard et al. 
2003, Duchene et al. 2006 and references therein), so its proper motion is likely to be affected 
by the gravitational influence of the other members of the system. As a consequence, the 
motion is likely to be curved and accelerated, rather than linear and uniform. To take this 
into account, we have also made fits to the data that include acceleration terms. This leads 
to functions of the form: 

a(t) = ao + (n a o cos 5)t + -(a a cos 5)t 2 + tc f a (t) (5) 
5(t) = 5 + fisot +^a s t 2 + irf s {t), (6) 

where fi a0 and fiso are the proper motions at a reference epoch, and a a and as are the 
projections of the uniform acceleration. Note that the acceleration undergone by a body in 
Keplerian orbit is usually not uniform. Assuming a uniform acceleration is acceptable here, 
however, because our data cover only a small portion (~ 2 yr) of the orbital period (a few 
decades -Duchene et al. 2006) of T Tau Sb. If VLBA data are obtained regularly in the next 
few decades, a full orbital fit will become possible, and indeed, necessary. 

The astrometric parameters were determined by least-square fitting the data points 
with either Eqs. [TH2] or Eqs. [5H6] using a Singular Value Decomposition (SVD) scheme 
(see Appendix for details). To check our results, we also performed two other fits to the 
data, a linear one based on the associated normal equations, and a non-linear one based 
on the Levenberg-Marquardt algorithm. They gave results identical to those obtained us- 
ing the SVD method. The reference epoch was taken at the mean of our observations (JD 
2453233.586 = J2004.627). 

4. Results 

The fit to the data points by Eqs. [TH2] (Fig. la) yields the following astrometric param- 
eters: 



aj2004. 6 27 = 04 h 21 m 59?425081 ± 05000005 

£j2004. 6 27 = 19°32'05"71566 ± 0"00003 

fi a cos5 = 4.00 ± 0.12 mas yr -1 

[is = —1.18 ± 0.05 mas yr -1 
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7T = 6.90 ± 0.09 mas. 



This corresponds to a distance of 145 ± 2 pc. The post-fit rms, however, is not very good 
(particularly in right ascension: ~ 0.2 mas) as the fit does not pass through many of the 
observed positions (Fig. la). As a matter of fact, 75 micro-arcsecond and -most notably- 
16.5 microseconds of time had to be added quadratically to the formal errors listed in Tab. 1 
to obtain a reduced x 2 °f 1 m both right ascension and declination; the errors on the fitted 
parameters quoted above include this systematic contribution. These large systematic errors 
most certainly reflect the fact mentioned earlier that the proper motion of T Tau Sb is not 
uniform because it belongs to a multiple system. Indeed, the fit where acceleration terms 
are included is significantly better (Fig. lb) with a post-fit rms of 60 yuas in right ascension 
and 90 //as in declination. It yields the following parameters: 



«J2004.627 
$J2004.627 
/V J2004.627 C °S 5 
I^S, J2004.627 

a a cos 5 
a s 

7Y 



04 h 21 m 59. s 425065 ± 0f000002 
19°32'05"71566 ± 0"0004 

4.02 ± 0.03 mas yr" 1 

-1.18 ± 0.05 mas yr~ x 

1.53 ± 0.13 mas yr~ 2 

0.00 ± 0.19 mas yr" 2 

6.82 ± 0.03 mas. 



To obtain a reduced x 2 °f 1 m both right ascension and declination, one must add quadrat- 
ically 3.8 microseconds of time and 75 microseconds of arc to the statistical errors listed in 
Tab. 1. The uncertainties reported above and in the rest of this article include this system- 
atic contribution. Note also that the reduced x 2 f° r t ne fit without acceleration terms is 
almost 8, if the latter systematic errors (rather than those mentioned earlier) are used. 

The trigonometric parallax obtained when acceleration terms are included, corresponds 
to a distance of 146.7 ± 0.6 pc, somewhat larger than, but consistent within 1.5<r with the 
value reported by Loinard et al. (2005). Recall, however, that this 2005 result was based on 
data that had been corrupted by a problem in the VLBA correlator; we consider the present 
value significantly more reliable. The present distance determination is somewhat smaller 
than, but within la of the distance obtained by Hipparchos (d = 177tgg pc). Note that the 
relative error of our distance is about 0.4%, against nearly 30% for the Hipparchos result, a 
gain of almost two orders of magnitude. 
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5. Implications for the properties of the stars 

Having obtained an improved distance estimate to the T Tauri system, we are now in a 
position to refine the determination of the intrinsic properties of each of the components of 
that system. Since the orbital motion between T Tau N and T Tau S is not yet known to 
very good precision, we will use synthetic spectra fitting to obtain the properties of T Tau 
N. For the very obscured T Tau S companion, on the other hand, we will refine the mass 
determinations based on the orbital fit obtained by Duchene et al. (2006). 



5.1. T Tau N 

The stellar parameters (T e g and -^ do j) of T Tau N were obtained by fitting synthetic 
spectra (Lejeune et al. 1997) to the optical part of the spectral energy distribution. In the 
absence of recently published optical spectra with absolute flux calibration, we decided to 
use narrow-band photometry taken at six different epochs from 1965 to 1970 (Kuhi 1974). In 
order to eliminate the contamination by the UV/blue (magnetospheric accretion) and red/IR 
(circumstellar disk) excesses, we restricted the fit to the range 0.41-0.65 /im. Two points 
at AA 0.4340, 0.4861 /im display large variations between epochs, they were also discarded 
as they are likely to be contaminated by emission lines. As a consequence, 56 photometric 
measurements at 13 wavelengths and 6 epochs had to be fitted. (See Fig. 2). We assumed that 
the star kept constant intrinsic parameters over the 5 years of observation, but allowed the 
circumstellar extinction to vary. Such an hypothesis is supported by long-term photometric 
observations (1986-2003) that show color-magnitude diagrams of T Tau elongated along 
the extinction direction (Melnikov & Grankin 2005); Kuhi (1974) also measured significant 
extinction variation in the period 1965-1970 using color excesses. 



0.4 



0.6 0.8 



X (jim) 

Fig. 2. — Fit to the photometry at six different epochs. Black squares are point that have 
been fitted; white squares represent other wavelengths excluded from the fit. 
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The non-linear fitting procedure used the Levenberg-Marquardt method and the deter- 
mination of errors was done using a Monte Carlo simulation. The synthetic spectra were 
transformed into narrow-band photometry by integration over the bandwidth of the mea- 
surements (typically 0.05 pm). As the fitting procedure could not constrain the metallicity, 
we assumed a solar one. Several fits using randomly chosen initial guesses for T g g, L-^ ^, 
and extinctions were performed in order to ensure that a global minimum x 2 was indeed 
reached. The errors reported by Kuhi (1974; 1.2%) had to be renormalised to 5.9% in order 
to achieve a reduced x 2 °f 1- This could result from an underestimation by the author or 
from positive and negative contamination by spectral lines -indeed, Gahm (1970) reports 
contamination as high as 20% for RW Aur. The best least-squares fit is represented in Fig. 
3, and yields T e g = 5112+g7 K and L-^ Q \ = 5.11^^1 L Q . The extinction varies between 
1.02 and 1.34, within 1-a of the values determined by Kuhi (1974) from color excesses. The 
effective temperature is consistent with a Kl star as reported by Kuhi (1974). 

In order to derive the age and mass of T Tau N, pre-main-sequence isochrones by 
D'Antona & Mazzitelli (1996) and Siess et al. (2000) were used. The fitting procedure was 
identical to the previous one: the age and mass were converted into effective temperature and 
luminosity, which in turn were converted into narrow-band photometry using the synthetic 
spectra. The derived parameters are shown in Tab. 3. The masses (1.831°;^ and 2.14TJq;}o 
M ) have overlapping error bars and are consistent with values found in the literature (e.g 
Duchene et al. 2006). The predicted ages, on the other hand, differ by a factor of 2. While 
the isochrones by D'Antona & Mazzitelli (1996) give an age in the commonly accepted range 
(1.15±g:ii M y r )' a somewhat larger value (2.39^;^) is derived from Siess et al. (2000). Note 
that the errors on the derived parameters are entirely dominated by the modeling errors; the 
uncertainty on the distance now represents a very small fraction of the error budget. 
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Aa [mas] Aa [mas] 

Fig. 3. — (a) VLB A positions (red squares) registered to T Tau Sa overimposed on the 
elliptical fit proposed by Duchene et al. (2006). Also shown are the velocity and acceleration 
vectors for our mean epoch deduced from our observations, as well as their counterparts 
from the fit by Duchene et al. (shown in blue). The dotted black lines around the measured 
acceleration and velocity show the error cone on the direction of each of these vectors, (b) 
Zoom on the region corresponding to our observations. In addition to the orbit and the 
VLBA positions, we show (in red) our best parabolic fit to our positions, as well as several 
recent infrared positions (in blue). The 2003.881 position is from Duchene et al. (2005), the 
2004.967 and 2005.868 positions are from Duchene et al. (2006), and the 2005.930 position 
is from Schaefer et al. (2006). 
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5.2. T Tau S 

The two members of the T Tau S system have been studied in detail by Duchene and 
coworkers in a series of recent articles (Duchene et al. 2002, 2005, 2006). The most massive 
member of the system (T Tau Sa) belongs to the mysterious class of "infrared companions" , 
and is presumably the precursor of an intermediate-mass star. T Tau Sb, on the other hand 
is a very obscured, but otherwise normal, pre-main sequence Ml star. The mass of both T 
Tau Sa and T Tau Sb were estimated by Duchene et al. (2006) using a fit to their orbital 
paths. Those authors used the distance to T Tauri deduced from Loinard et al. (2005). 
Using the new distance determination obtained here, we can re-normalize those masses. We 
obtain M Sa = 3.10 ± 0.34 M , and M Sb = 0.69 ± 0.18 M . These values may need to 
be adjusted somewhat, however, as the fit to the orbital path of the T Tau Sa/T Tau Sb 
system is improved (see below). Note finally, that the main sources of errors on the masses 
are related to the orbital motion modeling rather than to the uncertainties of the distance. 

6. Implications for the orbital motions 

T Tau Sb is a member of a multiple system, so it would be desirable to give its po- 
sition and express its motion relative to the other members of the system, T Tau N and 
-particularly- T Tau Sa. Since only T Tau Sb is detected in our VLBA observations, how- 
ever, registering the positions reported here to the other members of the system involves a 
number of steps. The absolute position and proper motion of T Tau N has been measured to 
great precision using over 20 years of VLA observations (Loinard et al. 2003), so registering 
the position and motion of T Tau Sb relative to T Tau N is fairly straightforward. Combin- 
ing the data used by Loinard et al. (2003) with several more recent VLA observations, we 
obtained the following absolute position (at epoch J2000.0) and proper motion for T Tau N: 

04 h 21 m 59. s 4321 ± OfOOOl 
19°32'06"419 ± 0"002 
12.35 ± 0.04 mas yr" 1 
-12.80 ± 0.06 mas yr _1 . 

Subtracting these values from the absolute positions and proper motion of T Tau Sb, we can 
obtain the positional offset between T Tau Sb and T Tau N, as well as their relative proper 
motion. For the median epoch of our observations, we obtain: 



«J2000.0 — 

^2000.0 — 

fi a cos S = 

M«5 = 
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/i Q cos5(Sb/N) = -8.33 ± 0.07 mas yr" 1 
/i 5 (Sb/N) = +11.62 ± 0.11 mas yr" 1 . 

The second step consists in registering the position and motion of T Tau Sb to the center 
of mass of T Tau S using the parabolic fits provided by Duchene et al. (2006). Here, both 
the proper motion and the acceleration must be taken into account. For the mean epoch of 
our observations, we obtain: 



// Q cos5(Sb/CM) = +0.3 ± 0.9 mas yr- 1 

^(Sb/CM) = +9.3 ± 0.8 mas yr" 1 

a Q cos5(Sb/CM) = +1.4 ± 0.2 mas yr" 2 

a 5 (Sb/CM) = -0.1 ± 0.3 mas yr" 2 . 



The last correction to be made is the registration of the positions, proper motions, and 
accelerations to T Tau Sa rather than to the center of mass of T Tau S. This is obtained by 
simply multiplying the values above by the ratio of the total mass of T Tau S (i.e. Ms a +Msb) 
to the mass of Sa. Using the masses given by Duchene et al. 2006), we obtain: 



fi a cos <5(Sb/Sa) = +0.4 ± 1.1 mas yr 1 

/^(Sb/Sa) = +11.4 ± 1.0 mas yr" 1 

a a cos 5(Sb/Sa) = +1.7 ± 0.2 mas yr -2 

a 5 (Sb/Sa) = -0.1 ± 0.3 mas yr" 2 . 

These two vectors are shown in Fig. 3 together with the VLBA positions registered to T 
Tau Sa, several recent infrared observations and the elliptical fit obtained by Duchene et al. 
(2006). The final error on the VLBA positions is the combination of the original uncertainty 
on their measured absolute position, and of the errors made at each of the steps described 
above. The final uncertainty is about 3 mas in both right ascension and declination, and is 
shown near the bottom right corner of Fig. 3b. 

Given the uncertainties, the position of the VLBA source is generally in good agreement 
with the infrared source position measured at similar epochs. Indeed, the first 2 VLBA 
observations were obtained almost exactly at the same time as the infrared image published 



16 



by Duchene et al. (2005), and the positions match exactly. The position of the VLBA source 
at the end of 2004 is also in agreement within la with the position of the infrared source 
at the same epoch reported by Duchene et al. (2006). The situation at the end of 2005, 
however, is somewhat less clear. Extrapolating from the last VLBA observation (~ 2005.5) 
to the end of 2005 gives a location that would be in reasonable agreement with the position 
given by Schaefer et al. (2006) but clearly not with the position obtained by Duchene et al. 
(2006). Note, indeed, that the two infrared positions are only very marginally consistent 
with one another. 

Our VLBA observations suggest that T Tau Sb passed at the westernmost point of its 
orbit around 2005.0, whereas according to the fit proposed by Duchene et al. (2006), this 
westernmost position was reached slightly before 2004.0. As a consequence, the trajectory 
described by the VLBA source is on average almost exactly north-south, whereas according 
to the fit proposed by Duchene et al. (2006), T Tau Sb is already moving back toward the 
east (Fig. 3). We note, however, that the fit proposed by Duchene et al. (2006, which gives 
an orbital period of 21.7 ± 0.9 yr) is very strongly constrained by their 2005.9 observation. 
Shaefer et al. (2006), who measured a position at the end of 2005 somewhat more to the 
north (in better agreement with our VLBA positions), argue that they cannot discriminate 
between orbital periods of 20, 30 or 40 yr. Orbits with longer periods bend back toward the 
east somewhat later (see Fig. 10 in Schaefer et al. 2006), and would be in better agreement 
with our VLBA positions. 

Another element that favors a somewhat longer orbital period is the acceleration mea- 
sured here. According to the fit proposed by Duchene et al. (2006), the expected transverse 
proper motion and acceleration are (G. Duchene, private communication): 



Thus, while the expected and observed proper motions are in good agreement, the expected 
acceleration is significantly larger that the observed value (see also Fig. 3). A smaller value 
of the acceleration would be consistent with a somewhat longer orbital period. 

In summary, our observations appear to be in reasonable agreement with all the pub- 
lished infrared positions obtained over the last few years, except for the 2005.9 observation 
reported by Duchene et al. (2006). As a consequence, our data favor an orbital period 



fi a cos<5(Sb/Sa) 

MSb/Sa) 
a a cos5(Sb/Sa) 
a^Sb/Sa) 



= +1.7 ± 0.2 mas yr" 1 

= +12.1 ± 1.2 mas yr~ 

= +3.1 ± 0.5 mas yr' 2 

= -1.4 ± 0.2 mas yr~ 2 



-l 
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somewhat longer than that obtained by Duchene et al. (2006). Exactly how much longer is 
difficult to assess for the following reason. The orbit proposed by Duchene et al. (2006) was 
obtained by fitting simultaneously the observed positions with a superposition of an elliptical 
path (of Sb around Sa) and a parabolic trajectory (of Sa around N). As a consequence, a 
modification of the Sa/Sb elliptical orbit (as may be required by our data) will result in a 
change in the parameters of the parabolic fit. But we use the latter to register our VLB A 
positions, proper motions, and accelerations against T Tau Sa. Thus, an entirely new fit will 
be needed to take into account the present VLBA observations. Such a fit will be presented 
in a forthcoming paper, where the numerous VLA observations available from the literature, 
as well as the VLBI observation from Smith et al. (2003) will also be taken into account. 



7. Conclusions and perspectives 

Using a series of 12 radio-continuum VLBA observations of T Tau Sb obtained roughly 
every two months between September 2003 and July 2005, we have measured the trigono- 
metric parallax and characterized the proper motion of this member of the T Tauri multiple 
system with unprecedented accuracy. The distance to T Tau Sb was found to be 146.7 ± 0.6 
pc, somewhat larger than the canonical value of 140 pc traditionally used. Using this precise 
estimate, we have recalculated the basic parameters of all three members of the system. The 
VLBA positions are in good agreement with recent infrared positions, but our data seem to 
favor a somewhat longer orbital period than that recently reported by Duchene et al. (2006) 
for the T Tau Sa/T Tau Sb system. 

Finally, it should be pointed out that if observations similar to those presented here 
were obtained regularly in the coming 5 to 10 years, they would greatly help to constrain 
the orbital path (and, therefore, the mass) of the T Tau Sa/T Tau Sb system. 

L.L., R.M.T, L.F.R., and R.A.G. acknowledge the financial support of DGAPA, UNAM 
and CONACyT, Mexico. NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. We are indebted to Gaspard 
Duchene for calculating the expected velocity and acceleration from his fit, and for his 
comments on the manuscript. We also thank the anonymous referee for his/her constructive 
comments on this paper. 



A. Appendix 



The parameters determined in this article (position at a reference epoch, trigonometric 
parallax, proper motions and accelerations) were obtained by minimizing the sum (Xa + Xs) 
of the residuals in right ascension and declination. The corresponding general mathematical 
problem is that where two functions x and y depend linearly on N independent parameters 
(dj and bi for x and y, respectively), and M common parameters cf. 



N M 
X(t) = ^OiUifr) +^2cjW*(t) 
i=l j=l 
N M 

y(t) = J2hv t (t)+Y,c J w y 1 (t) 

i=l j=l 

The values Xk and y k of the functions x and y have been measured at P times t k with errors 
<j£ and a y k) respectively, and the total \ 2 can De written: 



X = 



2 , 2 
Xx ' Xy 



P 

E 

fc=i 



-, 2 



+ 



ol 



Defining the following matrix elements: 



Ui(t k ) 



Vi(t k ) 



Wjitk) 



a 



x k 
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, Vk 
Wk = — , 

°k 



we can re-write the total x given by equation IA1I as 



X 



E 

k=l 



N M 

Ok - OtkiCLi - 22 1 k 3 C 3 

i=l 3=1 



+ 



AT M 
3=1 



i=l 



This sum of quadratic terms can clearly be seen as the squared norm of the vector: 



/ 01 - Sill a U a i ~ E/=l Jlj C 3 \ 
°2 - Ei=l "2ifli - Ej=l l2jCj 

4>i ~ Eii - Ei=i 5yCj 



Re-arranging the terms in this expression, we can re-write the total \ 2 as: 



X 



a 7 
0/3 5 







(A2) 



The procedure then consists of finding the vector X that minimizes an expression of the 
form: 

WA.X-BW 2 



An efficient algorithm to perform this operation is known as the Singular Value Decomposi- 
tion (SVD -see Press et al. 1992). This method is based on the linear algebra theorem that 
states that any I x J rectangular matrix M whose number / of rows is larger or equal to its 
number J of columns, can written as the product of an / x J column-orthogonal matrix U 
by a J x J diagonal matrix W with positive or zero elements, and the transpose of a J x J 
orthogonal matrix V: 



M = U W V T 



(A3) 



The rectangular matrix in expression IA2I has 2P rows (=24 in our case) and 2N + M 
column (=5 or 7, for uniform or accelerated proper motions, respectively) and can clearly be 
decomposed in that fashion. Since both U and V in the previous expression are orthogonal, 
their inverses are just their transposes. Also, if none of its diagonal elements are zero (which 
will be the case in all situations considered here), the inverse of W is a diagonal matrix 
whose elements are just the inverses of those of W. Thus, the inverse of matrix M can be 
written as: 

M~ l = V W~ l U T 

It can be shown (see Press et al. 1992) that, if the matrix M can be decomposed as above, 
then the vector X that minimizes the expression ||AX — B\\ 2 is simply: 

X = V W' 1 U T B (A4) 

An efficient way of solving our least-squares fit problem is, therefore, to form the rectangular 
matrix that appears in expression IA21 decompose it as in \A3>\ and calculate the value of the 
cij's, bi's, and q's using lA~4l This method was implemented in FORTRAN following Press et 
al. (1992). 

REFERENCES 

Andre, P., Deeney, B. D., Phillips, R. B., & Lestrade, J.-F. 1992, ApJ, 401, 667 

Beasley, A. J., Gordon, D., Peck, A. B., Petrov, L., MacMillan, D. S., Fomalont, E. B., & 
Ma, C. 2002, ApJS, 141, 13 

Bertout, C, Robichon, N., & Arenou, F., 1999, A&A, 352, 574 

Condon, J. J. 1997, PASP, 109, 166 

D'Antona, F., & Mazzitelli, I. 1996, ApJ, 456, 329 

Duchene, G., Beust, H., Adjali, F., Konopacky, Q. M., & Ghez, A. M. 2006, A&A, 457, L9 
Duchene, G., Ghez, A.M., & McCabe, C, 2002, ApJ, 568, 771 



- 21 - 



Duchene, G., Ghez, A. M., McCabe, C, & Ceccarelli, C. 2005, ApJ, 628, 832 

Elias, J. H. 1978, ApJ, 224, 857 

Elias, J. H. 1978, ApJ, 224, 453 

Evans, N. J., II, et al. 2003, PASP, 115, 965 

Feigelson, E.D., & Montmerle, T., 1999, ARAA, 37, 363 

Fomalont, E. B. 1999, Synthesis Imaging in Radio Astronomy II, 180, 463 

Gahm, G. F. 1970, ApJ, 160, 1117 

Giidel, M., et al. 2007, A&A, 468, 353 

Johnston, K.J., Gaume, R.A., Fey, A.L., de Vegt, G, & Claussen, M.J., 2003, AJ, 125, 858 
Kenyon, S.J., Dobrzycka, D., & Hartmann L., 1994, AJ, 108, 1872 
Knude, J., & Hog, E. 1998, A&A, 338, 897 
Kuhi, L. V. 1974, A&AS, 15, 47 

Lejeune, T., Cuisinier, F., & Buser, R. 1997, A&AS, 125, 229 
Mel'Nikov, S. Y., & Grankin, K. N. 2005, Astronomy Letters, 31, 427 

Loinard, L., Mioduszewski, A. J., Rodriguez, L. F., Gonzalez, R. A., Rodriguez, M. I., & 
Torres, R. M. 2005, ApJ, 619, L179 

Loinard, L., Rodriguez, L. F., D'Alessio, P., Rodriguez, M. I., & Gonzalez, R. F. 2007, ApJ, 
657, 916 

Loinard L., Rodriguez, L.F., & Rodriguez, M.I., 2003, ApJ, 587, L47 
Perryman, M.A.G, Lindegren, L., Kovalevsky, J., et al., 1997, A&A, 323, L49 
Phillips, R. B., Lonsdale, C. J., & Feigelson, E. D. 1993, ApJ, 403, L43 
Pradel, N., Chariot, P., & Lestrade, J.-F. 2006, A&A, 452, 1099 

Press, W. FL, Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes, 
Cambridge: University Press, 1992, 2nd ed. 

Schaefer, G. H., Simon, M., Beck, T. L., Nelan, E., & Prato, L. 2006, AJ, 132, 2618 



-22 - 



Seidelmann, P. K., & Fukushima, T. 1992, A&A, 265, 833 
Siess, L., Dufour, E., & Forestini, M., A&A, 358, 593 
Skinner, S. L., & Brown, A. 1994, A J, 107, 1461 

Smith, K., Pestalozzi, M., Giidel, M., Conway, J., & Benz, A.O., 2003, A&A, 406, 957 

Thompson, A. R., Moran, J. M., & Swenson, G. W. 1986, Synthesis Imaging in Radioas- 
tronomy, New York, Wiley-Interscience, 1986 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Table 2. Julian dates and Earth coordinates 



Mean UT date 


JD 


Earth Barycentric coordinates 


(yyyy.mm.dd hh.mm) 
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Table 3: Parameters of T Tauri N 

Isochrone set Siess et al. (2000) D'Antona & Mazzitelli (1997) 
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